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ABSTRACT 

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, 
cost, and risk of space-bome optical systems in the presence of dynamic disturbances. More robust systems 
almost certainly will be a result as well 

Many future space-based and ground-based optical systems will employ optical control systems to enhance 
imaging performance. The goal of the optical control subsystem is to determine the wave ont a errations 
and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently 
clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These 
elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable 
mirrorfs), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic 
control. These control subsystems may become elaborate indeed. But robust system performance will 
require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate 
systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or 

atmospheric sensors for dynamic imaging of individual, severe storms. The techno ogy eve ope ere 

could have a substantial impact on the development of new systems in the Space Science Enterprise, sue 
as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric 
Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic 
Structure, could benefit. These systems most likely will contain large, flexible optomechanical structures 
subject to dynamic disturbance. Furthermore, large systems for high resolution imaging o p anets or e 
sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small 
targets as well and could benefit from the technology developed here. 

We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static 
aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should 
work for all but the lowest light levels and highest frequency dynamic environments. This technique may 
serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk 
operation. For a given object, a short exposure image is •‘frozen’’ on the focal plane in the presence of the 
environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to- 

frame linear shift invariance. Therefore, although the Point Spread Function is vary ing from frame to 

frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limi 
of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. he 
magnitude of the complex object spectrum is straightforward to determine by well-established approaches. 
The phase of the complex object spectrum is decomposed into two pans. One is a single-valued function 
determined bv the divergence of the optical phase gradient. The other is a multi-valued function determine 
bv the circulation of the optical phase gradient-"hidden phase.” Finite difference equations are developed 
for the phase. The noveltv of this approach is captured in the inclusion of this -hidden phase. his 
technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure 
frames while simultaneouslv estimating the phase as a function of time from a set ot exposures. 



The following subsections of the paper are included: Introduction. Technical Approach. Plans. Preliminary 
Results of earlier work and a Conclusion. 


1. Introduction 

The desired result of this new algorithm is the development and verification in a laboratory setting of a 
technique to reduce the sensitivity and improve performance of NASA and DOD optical systems in 
dynamic environments. On the ground and in the laboratory such systems are sensitive to an t e 
performance is degraded by atmospheric turbulence, laboratory seeing, and motion/jitter ot key optica 
elements. In space the sensitivity of such optical systems is primarily to jitter; however some effects o 
downward looking through an atmosphere may still be presence. Of course, seeing effects may still be 
present if the system is sealed and operating at non-zero pressure. In any event the development envisione 
here will be able to improve algorithms that compute system performance. A key feature ot the 
improvement lies in the inclusion of the “hidden phase” component in this algorithm. 

We intend to investigate the improvement, our technique offers, in the context of various system 
architectures. For large, space systems with flexible structures various control architectures wi l be 
applied. Control will most likely be distributed over several elements. These elements can mclud 
primary mirror segments, secondary, tertiary, deformable mirror(s)-DM(s), Fast SteemgMirTorisKS^), 
etc. The control subsystem goal is to improve a system with substantially aberrated l imaging 
limited imaging. Further the control subsystem will have to ameliorate the effects of temporally an 
sZSyTa^g disturbances. I. is in .he reahn of the tempontlly vtnying disnnbances that thts algomta 
will be most effective. Indeed it is the hope that it can eliminate the need for any higher order aberration 
correction in the temporal domain and leave perhaps just a simple two-axis FSM. Post-processmg usmg 
this algorithm will work at all but the lowest light levels. Realistic estimates of light-level limits on 
performance will be investigated. 

Space systems of the future will be driven to lighter weight and lower areal de ™ d ^ 

systems tend to be flexible and therefore prone to dynamic/jitter effects. For all NASA and y , 

sources of dynamical disturbances/jitter are certainly reaction wheels (as well as other stabilization devices), 

detector coolers, as well as other sources. 

The Severe Storm type of mission may be particularly exciting and challenging, since the i object of study 
(a single storm cell ) is time evolving. The nature of the temporal disturbances onboard the satellite and 
those in the object will have to be studied and separated. Other geostationary systems may also benefit, 
such fourier transform spectrometers and other advanced meteorological satellites. Low earth or i g 
systems such as Landsat may benefit as well. 

The technology developed here could have a substantial impact on the development of new space 
telescopes- such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large 
Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter of *e Ev^t.on 

of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanica 
structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging o 
planets and the sun would like benefit from this technology. 

Positive results from this approach could lead directly to a high accuracy (near diffraction-limited) system 
for imaging in the presence of various disturbances, including vibration. Further our approach has the 
potential to reduce cost and risk for large, light-weight systems in many potential applications- 
approach is innovative and requires development to enhance its technology readiness. But it is believed 

hold great promise for improvements. 



For high power level systems, mainly DOD tactical and strategic systems, coolant flow induced jitter may 
be a significant source of disturbance. We have enumerated some specific types of systems that may be 
candidates for improvement. 

Flowing from this development may be significant reduction in the complexity of space-borne optical 
systems. Consequent significant reductions in cost may also accrue form such a reduction in complexity. 
More robust systems may be a result as well. The output from this development quite likely would be a 
quantitative and qualitative reduction in risk for such space missions. 


2. Technical Approach 

This section of the proposal is divided into three sections. In the first section we present an overview of the 
speckle imaging technique to be used in this investigation. Secondly we elaborate the specifics of the Knox- 
Thompson approach used here and explain the enhanced feature; the development of the hidden phase. In 
the last section we discuss the particular plan to be followed in the future. 


2.1 Overview 

Speckle imaging techniques have been evolving since the fundamental idea was presented almost 30 years 
ago.' Labeyrie’s key observation was that the speckles in short exposure images contain richer spatial ^ 
frequency information than long exposures. The speckle imaging technique can be described as follows. 
Two data sets are required to perform this imaging; one set of RO short exposures images of the object 
under study and one set RR of similar images of a bright, nearby reference. Normally, values of RO and 
RR are on the order of tens or hundreds, driven by signal-to-noise ratio considerations. Exposure times are 
typically a few to tens of milliseconds. 

The first step in the process is to compute the fourier transform of the images. Next one computes the 
modulus squared and the cross spectrum or bispectrum. We use the cross spectrum obtained by the Knox- 
Thompson method. Statistics are then accumulated and the process is repeated either RO or RR times as 
appropriate. Then one computes the average modulus squared and the average cross spectrum or 
bispectrum. Continuing, one deconvolves the image modulus squared and computes the phase spectrum. 
Finally the fourier modulus and phase data are combined; the inverse transform is computed; and the image 
estimate is obtained. 

Many methods exist to find the modulus and we will not elaborate them here. We will explain our approach 
below. Iterative methods start with the modulus and proceed along the following lines. First the inverse 
FFT is computed. Next the constraints are imposed. Finally the FFT is computed to obtain the phase and 
the process is iterated. We use a method to directly compute the phase using difference equations. 


2.2 Detailed Approach 

Here, we will concentrate on the particular implementation called the Knox-Thompson, or cross-spectrum, 
method and, in particular, deal with the object spectrum phase reconstruction. 


We label the arrays of short-exposure data as: 


The index 1 = 1.2, ... R labels the frame, while the next two indices, (ij). label the pixel in the array. We 
take a two-dimensional Finite-Fourier-Transform (FFT) on the pixel indices as. 



( 2 ) 


d{i,m,n)=FFT{d(l\ij)}. 

The array indices, m and n, label the elements of the two-dimensional FFT. The zero spatial frequency 
point is labeled ( m , n ) . We then form the unbiased cross-spectra as 

C t (w,«) = -^y^ | d{l,m,n)d (< l\m + \,n)-d (/;m+l,«)} 


C y {m,ri) = ^^ \d{i,m,n)d*(l;m,n + \)- d (l\m,n+ l)J. 

/=1 

in which we have picked one unit offset in each direction in the FFT; other offsets can easily be analyzed, 
however, the unit offset choice gives the best accuracy on the high-spatial-frequency phase. The ensemble 
averaging over many realizations of the turbulence-induced phase screens makes the cross-spectrum trans er 
functions real-valued. 3 Therefore, the turbulence-induced phase errors are averaged away and the phase of 
the cross-spectrum can be directly related to phase differences in the object spectrum. The remaining 
phases of these two cross-spectra lead to two two-dimensional difference equations for the phase of the 
object spectrum as: 


C x {m, n) = \C x {m, n) | exp( i V x ) 


C y (m,n)=\Cy{m,n) expf/Fy} 


(4a) 


in which the cross-spectrum phases, V x , V y , are related to the object spectrum phase, W, as. 


V x (m,n) = W(m,n)-W(m + l,n) 
Vy (m, n)=W (m, n) - W (m , n + l). 


(4b) 


In a new paper 4 these equations for the object spectrum phase are solved 4 ; we will refer to these equations 
as "phase flow” equations. In the derivation 4 we first define discrete gradient, divergence and curl 
operations, followed by our proof that the phase flow on the FFT-grid can be decomposed into an 
irrotational part with zero curl and a rotational part with zero divergence. Next in the derivation we wi 
demonstrate that in a similar fashion, the object spectrum phase can be decomposed into a regular single- 
valued function determined by the divergence of the phase gradient, as well as a multi-valued function 
determined by the circulation of the phase gradient; this second function has been called the “hidden 
phase”. 5 We will present a solution method that gives both the regular and hidden parts of the object 
spectrum phase. Elsewhere it is also demonstrated 5 that the standard least-squares solution to the two- 
dimensional difference equations will only generate the regular part of the phase, while always missing the 
hidden part. In the preliminary results section , we will give several examples of imaging through 
turbulence and post-processing, augmenting the Knox-Thompson method with a hidden phase algorithm. 

In particular, we will compare reconstructions with and without the hidden phase component included in the 
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reconstruction algorithm. Finally, in the conclusion we discuss prospects for the method, as well as ongoing 
problem areas. 

The solution difficulties associated with the cross-spectrum phases are reasonably well-known. In 1986^ 
Fontanella, et. al. provided a discussion of the effects of phase dislocations on the phase reconstruction. In 
1987. Takajo, et. al. acknowledged the difficulty of phase reconstruction near isolated zeroes in the object 
spectrum and proposed a solution. 7 More recently , several authors have presented branch-point to erant 
phase reconstructors. 8 - 910 The present paper 4 adds to these results by emphasizing a unique decomposition 
of the phase flow into an irrotational and a rotational part. We feel that this approach lends considera e 
claritv to the phase reconstruction problem. Furthermore, we specifically treat phase reconstruction in t e 
Knox-Thompson speckle imaging method, providing a post- processing algorithm that appears to be very 
effective for ground-based astronomy applications. 


We have mentioned above that we can decompose the vector phase flow into a piece with zero divergence, 
and a piece with zero circulation . That a continuous vector field can be decomposed into rotational an 
irrotational parts has proved quite useful in hydrodynamics for nearly 150 years. Our results are a discrete 
analog of those well-known results. 


2.3 Plans 


We now present the future direction for conducting this development. It can be decomposed into four 
phases. These phases are enumerated below: 

(1) Complete application of the method to the atmospheric turbulence case 

(2) Develop disturbance model for the space environment 

(3) Obtain calculated results for this environment 

(4) Develop designs and plans for laboratory demonstrations of the technique. 

The first phase will entail expanding on results already obtained for atmospheric compensation and to be 
briefly described below. This work will include assessing precision of the method and developing light 
level/dynamic range models. Other source objects including distributed ones will be investigated. The 
modeling of laboratory (horizontal) seeing effects will commence 


The second phase will cover the development of the jitter model for two candidate prototype systems; one a 
two element telescope system and the other a large observatory with three or four mam telescope elements. 
Jitter sources on spacebome systems will be investigated such as gyros, stabilization subsystems, cryogenic 
subsystems, cooling subsystems, etc. Disturbance sensitivities of the optical subsystem will be explored, 
such as motion of key optical elements (secondary, tertiary, etc.) and the detector. The impact of the 
statistical distribution functions (Gaussian, Poisson, etc.) will be incorporated. Lab seeing effects will e 
fully incorporated. Development of typical system transfer models will start. In these models we wi 
determine the attenuation of the disturbances expected with candidate control systems and the resulting 
residual after control. 


In the third phase we will fully exercise the baseline algorithm on typical source objects with known jitter 
disturbances present. Image reconstruction performance of the algorithm will be thoroughly evaluated. 
A^ain we will assess precision of the method and develop light leveL'dynamic range models. Other source 
objects including distributed ones will be investigated as before for atmospheric compensation. System 
transfer performance will be fully evaluated for the two specific candidate systems mentioned in the 
paragraph above. Comparison of this speckle-imaging algorithm to other such approaches will be starte . 
The relationship of speckle approaches to phase retrieval and phase diversity reconstruction approaches wi 
be addressed. 


In the final phase, a laboratory demonstration concept will be developed. A simple, probably two 
element, telescope system will be designed and analyzed. A surrogate source will be image on a 
CCD detector. Stingers will be identified for the telescope secondary and the detector to introduce jitter. 



1 he performance of (he image reconstruction process in the presence ot known disturbances an 

Known sensitivities will be analyzed The comparison of this speckle imaging algorithm to other 

such approaches will be finished. The relationship of speckle approaches to phase retrieval and 

phase diversity reconstruction approaches will be further addressed. Avenues to obtain actual hardware and 

lab space will be pursued. 


3. Preliminary Results 

In this section, we will demonstrate image reconstructions, using the Knox-Thompson Method augmented 

bv our complete solution for both the regular and hidden parts of object spectrum phase. 

Figure (la) plots a diffraction-limited image of a four-star group, in which one ot the stars is twice : e 

intensitv of the other three. (We show all figures on a gray scale on which the original j2 x o intensity data 
is bilinearlv interpolated to a 256 x 256 grid. The specific bilinear interpolation method is in Ret. 1-.) 
Figure (lb) shows a turbulence-degraded, long-time exposure of the identical astensm. Here we l av 
averaged 640 short exposures, using aberrated point-spread functions. In addition we have added sh^ 
noise to each pixel (32 x 32 array) in each short-exposure frame, while assuming 10,000 photons per frame. 



We use standard speckle methods, as described in reference 1, to determine the modulus of the o^ectFFT. 
In this step we use an independent set of reference star short exposures, in order to deconvolve the speckle 
Iranstl Action from the modulus of the object FFT. We then use the Knox-Thompson method to generate 
difference equations for the phase, and to solve for both the regular. W . and multi-valued. W , 
components of the object FFT phase. 

Prior to taking the inverse FFT. we multiply by the optical transfer (unction. H DL (m.m. of the diffraction- 
limited telescope. The final diffraction-limited object reconstruction is then given as 


/(/.;)= F/T -1 j e 


W S {m,n)+W H (/».*)) H DL (m.n)\ . 


This reconstruction is plotted in Fig. ( 1 c); the reconstruction accuracy is excellent. It we affect the 
identical reconstruction procedure, but use a minimum least-squares phase estimator, then the h'ddenph^e 
is not included in the object FFT estimate. The reconstruction results with the hiduen phase neg e 
shown in Fie i Id). Obviously, the hidden phase plays a critical role in the reconstruction process. 


h 




Figure 1c.) Reconstruction with hidden phase. 


Figure Id.) Reconstruction without hidden phase. 


As a second example. Fig. (2a) shows a contour plot of a diffraction-limited image ot a binary star in whic 
one star is twice the intensity of the other. Figure (2b) shows the turbulence-degraded . long-time average 
of 640 short exposures with 10.000 photons per frame. Using the augmented Knox-Thompson procedure, 
we obtain the results in Fig. (2c) In this simpler object case, there are no zeroes in the original object 
spectrum, and the hidden phase is essentially zero. Figure (2d) shows the reconstruction when the hidden 
phase is set exactly to zero. The results are excellent in either case. 








Figure 2c.) Reconstruction with hidden phase. Figure 2d.) Reconstruction without hidden phase. 


4, Conclusion 


The Knox-Thompson, or cross-spectrum, method provides two two-dimensional difference equations for 
the phase of the object spectrum. A complete solution to these difference equations can be decomposed 
into a regular single-valued function, W s , determined by the divergence of the phase gradient, as well as a 
multi-valued function, W H , determined by the circulation of the phase gradient. Furthermore, this phase 
decomposition is unique and complete. An accurate object spectrum reconstruction then requires a modulus 
estimate, as well as both parts of the phase. For almost any object on the celestial sphere more complicated 
than an unequal magnitude binary, both components of the object spectrum phase are needed for an 
accurate reconstruction. Knox-Thompson reconstructions based on a minimum-least-squares phase 
estimator will usually be of limited value. 

The examples in Section 3 demonstrated reconstructions on 640 short-exposure frames in the presence of 
significant turbulence and shot noise. In addition, we have tested the method with realistic shot and read 
noise corrupting the short-exposure data. Although our augmented Knox-Thompson method worked well in 
these examples, as well as in many other test problems, several questions remain: Is the method an 
optimum method in the presence of noise? Is there a better way to motivate or derive our method tor 
filtering the phase flow circulation values? Can we quantify the reconstruction accuracy as a function of 
turbulence strength, photons per frame, and number of data frames 0 These questions, along with other 
evolving issues, will require further investigation. 


The technology developed here could have a substantial impact on the development ot new space 
telescopes; such as the Next Generation Space TelescopefiNGST) and its tollow-on the Next NGST. Large 
Interferometric S'. stems of non-zero field, such as Planet Finder and Submillimeter Probe ot the Evolution 
of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanical 
structures subject to dynamic disturbance. Furthermore, large systems tor high resolution imaging ot 
planets or Earth from space may also benefit. Such systems could include next generation Landsat systems 
or atmospheric sensors for dynamic imaging of individual, severe storms. Tactical and Strategic Defense 
sv stems will need to image very small targets as well and could benefit trom ihe technology developed here. 
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We believe our technique offers improvement in the context of various system architectures. For large, 
space systems with tlexible structures various control architectures will be applied. Control will most likely 
be distributed. The control subsystem goal is to improve a system with substantially aberrated imaging to 
diffraction-limited imaging. In the realm of the temporally varying disturbances our algorithm is expected 
to be the most effective. Indeed it is the hope that it can eliminate the need for any higher order aberration 
correction in the temporal domain and leave perhaps just a simple two-axis FSM. Post-processing using 
this algorithm will work at all but the lowest light levels. Realistic estimates of light-level limits on 
performance will be investigated. This algorithm could in the future be applied and the software written for 
a specific flight project. 


This investigation may yield significant reduction in complexity, cost, and risk (quantitative and qualitative) 
of space-borne optical systems by overcoming dynamic disturbances. More robust systems almost 
certainly will be a result as well. 

5, Acknowledgement 

We would like to acknowledge Michael K. Tiltons’s contributions in providing information and most 
particularly in developing the data and figures for this paper. 

6. References 

1. A. Labeyrie, “Attainment of diffraction limited resolution in large telescopes by Fourier analyzing 
speckle patterns in star images, " Astron. Astrophys. 6, 85-87 (1970). 

2. M. C. Roggemann and B. Walsh, Imaging Through Turbulence , CRC Press, 124-125, 1996. 

3. K. T. Knox and B. J. Thompson, “Recovery of images from atmospherically degraded short-exposure 
photographs," Astron. J. 193, L45-L48 (1974). 

4. G.C. Dente, “Speckle Imaging and Hidden Phase,” Appl. Opt. 39, 1480-1485 (2000). 

5. D. L. Fried, “Branch point problem in adaptive optics," J. Opt. Soc. Am. 15, 2759-2768 (1998). 

6. J. C. Fontanelia and A. Seve, “Reconstruction of turbulence-degraded images using the Knox- 
Thompson algorithm," J. Opt. Soc. Am. 4, 438-448 (1987). 

7. H. Takajo and T. Takahashi, “Least-squares phase estimation from the phase difference," J. Opt. Soc. 
Am. 5,416-425 (1988). 

8. W. W. Arrasmith, “Branch-point tolerant least-squares phase reconstructor," J. Opt. Soc. Am. A 16, 
1864-1872(1999). 

9. E. Le Bigot and W. J. Wild, “Theory of branch-point detection and its implementation," J. Opt. Soc. 
Am. A 16, 1724-1729 (1999). 

10. T. W. Lawrence, D. M. Goodman, E. M. Johansson and J. P. Fitch, “Speckle imaging of satellites ar 
the U. S. Air Force Maui Optical Station," Appl. Opt. 31, 6307-6321 (1992). 

11. A. Sommerfeld, “Mechanics of Deformable Bodies: Lectures on Theoretical Physics, Volume II , 
(Academic Press, New York, 1959), section IV. 20. 

12. W. H. Press. B. P. Flannery, S. A. Teukolsky and W. T. Vetterling, '‘Numerical Recipes in Fortran 77, 
Second Edition", (Cambridge University Press, 1992), page 1 16. 


9 



Richard G. Lyon, 10:45 AM 1/24/00 , Computational Optica Workshop 


; ° hn • - .Wo L c j.m j . L-i'jj s'-:. nas.i . 70 v, Howard . 
n woo L t 0 o;5 . a r;zo no . edu , ye Lena 0 ce sd L 3 L . g 3 
Sr ; w • '^ung . L igs fc . nasa . gov, vass ii is . z 
i< i ao L e 1 . Z ha pg. 1 0 gs f 0 . na 3 a . go 7 
Compu tat ior.a 1 Optics Workshop 


Wood . 10 gs t"c . rasa . gov, 

.nasa.gov, 

i f is0 Lmco . com. 




PLEASE FREELY DISTRIBUTE 


V+ilrlf + drV** 


/ 


■V * + + * ic k ★ + + 


+ + + ★ * f 


WORKSHOP ON COMPUTATIONAL OPTICS AND IMAGING FOR SPACE APPL 


CATIONS 


-ccstion: NASA/Goddard Space Flight Center 

Greenbelt, MD 20771 


Date 


Workshop - May 10, 11, 12, 2000 
500 Word Abstract /summary - March 24 
(please specify oral or poster presentation) 
4-8 page manuscript - May 10, 2000 


N-r 


Sponsored by NASA Goddard Space Flight Center, USRA F s Center of 
excellence in Space Data and Information Sciences and the Optical 
Society of America 


Scope of Workshop 

NASA Goddard Space Flight Center, USRA's Center of Excellence in Space 
Da^a and Information Sciences and the Optical Society of America are 
pleased to announce the first workshop on Computational Optics and Imaging 
for Space Applications. The general purpose of this meeting is to bring 
together researchers in all aspects of computational optics, imaging, 
on-board and ground processing and advanced computer architectures whose 
primary emphasis has been on space based applications. 

Background and Topics 

Space based imaging systems require increasingly larger aperture sizes to 
,<eep^pace with the demand for higher spatial and spectral resolution imagery 
p or ^arth and Space sciences missions. The cost and weight becomes 
increasingly prohibitive for telescopes and instruments with apertures greater 
than 1 meter. A number of solutions are possible and include - deployable 
segmented aperture systems, sparse aperture systems, interferometric imaging 
systems and computational deconvolution and super— resolution techniques. 

*ne^ commonality of these techniques lies in increased reliance on 
s v>pn_s _ i wa ;_ed computational and information theoretic techniques. We are 
nc^c^ng cnis workshop to address primarily the computational methods and 
-e--..r.*q^es required for each or the above and are soliciting both oral and 
pesters papers in the areas: 


Model ing/Simuiat ion of optical systems. 

- Wavefront Sensing and Optical Control. 

- Phase Retrieval and Phase Diversity. 

Image Restoration, Deconvolution and Super-Resol 
Sparse and Inter feromet r ic imaging algorithms. 

- computer aroni te ctures for or.-bcard and ground c 


*'i n _ _ e 


i s rea 


applications there is commons 1 
t : ;aoi^itate an exchanges of 
foster collaborative efforts. 


i c e a s 


a relatival 
overlap. Th 
and methods 



r .< s n ^ 


to peter. 


expect 
call v 


Printed for Eric Young <Enc . W . Young@gsf c . nasa . gov> 


2 


Rxchard G. Lyon, 10:45 AM 1/24/00 , Computational Optxca Workshop 


'h i i rman : 


Richard G. Lyon {NASA/ Go FC) 


Comm it tee : 


Dennis Andrucyk ( NA5A/GSFC) 
(NASA/GSFC) 
(NASA/GSFC) 
(ERIM) 
(NASA/GSFC) 
Dave Leisawitz (NASA/GSFC) 
Timothy Murphy (CESDIS/UMBC) 
John Wood (NASA/GSFC) Eric 


Chuck Bowers 
John Dorband 
James Feinup 
Bill Hayden 


Paul Atchescn (Bail) 

Scott Dixon (Photon Research Associates; 
Lee Feinberg (NASA/GSFC) 

‘ Milton Haiem (NASA/GSFC) 

Edward Howard (NOAA) 

Jacqueline Lemoigne (NASA/GSFC ; 

Richard Paxuian (ERIM) 

Young (NASA/GSFC) 


Contact Information 


Nora Farrell 
Georgia Flanagan 

Shelly Meyett 

Richard Lyon 


farrell@cesdis .gsfc.nasa.gov 301-286-1112 
georgia0cesdis . gsfc.nasa.gov 301-286-2080 

shell y@cesdis . gsfc. nasa . gov 301-286-8755 
lyon0jansky. gsfc. nasa . gov 301-236-4302 


A web page will shortly be setup for this workshop. 


Printed for Eric Young <Enc. W. Young^gsfc. nasa.gov> 


3 



